THE CORRELATION FUNCTION OF FLUX-LIMITED X-RAY CLUSTERS 



Neta A. Bahcall and Renyue Cen 
Princeton University Observatory 
Princeton, NJ 08544 USA 



To appear in 



The Astrophysical Journal Letters, May 1, 1994 



1 



ABSTRACT 

We show that the spatial correlation function of a flux-limited sample of X- 
ray selected clusters of galaxies will exhibit a correlation scale that is smaller than 
the correlation scale of a volume-limited, richness-limited sample of comparable 
apparent spatial density. The flux-limited sample contains clusters of different 
richnesses at different distances: poor groups are found nearby and rich clusters 
at greater distances. Since the cluster correlation strength is known to increase 
with richness, the flirx-limited sample averages over the correlations of poor and 
rich clusters. On the other hand, a volume-limited, richness-limited sample has a 
minimum richness threshold, and a constant mixture of richnesses with redshift. 
Using the observed correlation scale for rich (i? > 1) clusters, ro{R > 1) = 21 ± 
2h~^Mpc that was determined from previous volume-limited studies, we derive 
for the ROSAT flux-limited X-ray cluster sample Tq (flux— limited) ~ 14/i^^Mpc, 
in agreement with the recently observed value of 13.7 ± 2.3/i~^Mpc. 

Cosmology: large-scale structure of Universe - galaxies: clusters of - X-rays: 
galaxies 
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1. INTRODUCTION 

The correlation function of clusters of galaxies constrains cosmological mod- 
els for the formation and evolution of structure in the universe (White et al. 
1987; BahcaU 1988; Suginohara & Suto 1991; BahcaU & Cen 1992; Croft & Ef- 
stathiou 1993). The cluster correlation function is known to be stronger than 
the correlation function of galaxies; the correlation scale of the rich and rare 
clusters of richness class R > 1 is ro{R > I) = 21 ±2/j,~-'^Mpc (BahcaU & Soneria 
1983; Postman et al. 1992; Peacock & West 1992), while the galaxy correla- 
tion scale is ro{g) = 5h~^Mpc (Groth & Peebles 1977). The cluster correlation 
strength is observed to increase with cluster richness: rich, rare clusters exhibit 
stronger correlations than the more numerous poor clusters and groups (Bah- 
caU & Soneira 1983; BahcaU & West 1992). All samples of clusters studied so 
far, from poor groups of galaxies to the richest R > 1 and R > 2 Abell clus- 
ters, including the intermediate richness clusters observed by the new automated 
cluster surveys of the EDCC (Nichol et al. 1992) and APM (Dalton et al. 1992) 
are consistent with a universal richness-dependent correlation (BahcaU & West 
1992). Large-scale cosmological N-body simulations of galaxy clusters also show 
the same dependence of the cluster correlation function on richness (BahcaU & 
Cen 1992; Croft & Efstathiou 1993). This richness-dependent correlation func- 
tion applies to complete richness-limited samples of clusters (i.e., clusters above 
a given richness threshold); it represents the underlying spatial distribution of a 
system of clusters of a given richness class. 

Recently, the spatial correlation function of a sample of X-ray clusters of 
galaxies detected in a fliix-limited survey of the ROSAT X-ray satellite was re- 
ported (Romer et al. 1993). The sample contains all X-ray sources above a given 
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X-ray flux threshold that are associated with a local galaxy density enhance- 
ment. A correlation length of = 13.7 ± 2.3/i~^Mpc was determined for the 
sample. Romer et al.(1993) conclude, from a direct comparison of this correla- 
tion length and the larger length of the rich R> 1 clusters, that inconsistencies 
exist between the two results, and that the R > 1 cluster correlation scale has 
been overestimated. 

In the present letter we show that a flux-limited sample, such as the X-ray 
sample described above, differs significantly from a richness-limited sample (such 
as the R> 1 clusters from which the 21/i~^Mpc scale length was obtained). The 
flux-limited sample is a "richness-mixed" sample; it contains, by definition, poor 
clusters nearby and rich clusters at greater distances. Since cluster correlations 
depend on richness, a simple comparison between the correlation properties of a 
flux-limited and a richness-limited sample is inappropriate. Here we show that, 
based on the richness-dependent correlation function and a correlation scale of 
ToiR > 1) = 21 ± 2h~^Mpc for i? > 1 clusters, the expected correlation scale 
for the above flux- limited X-ray sample is To (fliix— limited) « 2/3ro{R > 1) = 
14:h~^Mpc, as is indeed observed. This flrst X-ray selected ROSAT cluster sample 
thus confirms the richness-dependent cluster correlation function. 

2. THE CORRELATION FUNCTION OF A FLUX-LIMITED 

X-RAY SAMPLE 

The richness-dependent cluster correlation function is represented by (Bah- 
caU & Soneira 1983; BahcaU & West 1992) 

ecc(r) ~ 4iVr-i-8 , (1) 

where ^cc = Ar~^-^ = {r/ro)~^'^ is the standard form of the correlation function 
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(with an amplitude A and a correlation scale Tq), and N is the median richness 
of the cluster sample. Relation (1) applies to complete volume- limited, richness- 
limited samples (all clusters above a threshold richness limit); it represents the 
spatial distribution of clusters of a given richness class. Similarly, the cluster 
correlation amplitude also depends on the mean separation of clusters, d (where 
d = n~^/^, and n is the space-density of the cluster sample). This led to the 
universal dimensionless cluster correlation function (Szalay & Schramm 1985; 
Bahcall & West 1992) 

^cc{d) = 0.2(r/d)-^-^ = (r/OAd)-^-^ , i.e., ^ OAd . (2) 

The above-described universal correlation function is seen both in observa- 
tions and in model simulations (Bahcall & Cen 1992). It applies to complete 
volume-limited, richness-limited samples, where n and d represent the underly- 
ing density and mean separation of a complete system of clusters above a given 
richness threshold. All the principal cluster samples analyzed to-date for which 
relations (1) and (2) apply are volume and richness limited (e.g., Abell, Zwicky, 
APM, EDCC clusters, and groups; see summary in Bahcall & West 1992). 

Recently, a flux- limited correlation function of X-ray clusters was determined 
by Romer et al. (1993). The sample includes all X-ray sources above a flux 
threshold of Fx > 10~^^ergs cm~^ s~^. A total of 161 sources associated with 
some enhancement in the galaxy density distribution are detected in a 3100 
deg^ region; 128 of these systems have measured redshifts. Romer et al. flnd 
To = 13.7 ± 2.3/?.~^Mpc for the correlation function of this flux-limited redshift 
sample (for a correlation slope of —1.9 ±0.4). They contrast this correlation scale 
with the value of = 21 ± 2/i~-'^Mpc observed for the rich R>1 clusters. 
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However, the volume- limited and the fliix-limited samples are not expected 
to yield the same correlation scales, since the criteria for inclusion in the two 
samples are different. The volume-limited sample includes all clusters above a 
given richness threshold; the flux-limited sample includes a mixture of richnesses 
that is a function of redshift. The average observed cluster density, n, and mean 
separation, d, of the flux-limited sample are not representative of a given richness 
system, and thus can not be applied in relations (1) and (2). 

What would the correlation function of such a flux-limited sample be if the 
underlying cluster correlation is represented by the richness-dependent universal 
correlation function (relations 1 and 2)? We address this question below. 

An observational relation exists, as expected theoretically, between the X- 
ray luminosity of clusters and cluster temperature, or mass. Henry & Arnaud 
(1991) find that, for a Hubble constant of Ho = 100/tkm s~^ Mpc~^ (where 
h = 1 wiU be used throughout), L^{Bo\) = 2.5 x io^'^t'^-7±oa^ ^^^^^ L^{Bo\) 
is the bolometric X-ray luminosity of the cluster in erg s~^, and T is the in- 
tracluster gas temperature in keV. Edge & Stewart (1991) find similar results, 
with L^(Bol) = 2.5 x lO^^T^-^ (see also Henry et al. 1992 and David et al. 1993 
for similar relations). Converting the bolometric luminosity to the ROSAT ob- 
served energy band of 0.1-2.4 keV (for the typical range of cluster temperatures 
T ~ 2 - 10 keV), we find La;(0.1 - 2.4keV) 2.8 x 10^2y2±o.4gj,gg rj.^^ 
virial mass of a cluster is proportional to the temperature, T, or to the square 
of peculiar velocity in the cluster, a: M oc T oc (Sarazin 1988; Bahcall & Cen 
1993; Lubin & BahcaU 1993). Lubin and BahcaU find, on average, a = 400T°-^ 
km s~^ for the best a oc T^'^ observed relation. The virial cluster mass within 
1.5/i~^Mpc of the cluster center, assuming an isothermal density profile, is then 
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M(< l.S/i-^Mpc) = 2a'^R{l.h)/G ^ 1.1 x 10^^T{keV). Combining the above 
relations we find 

L^{0.1 - 2.4keV) 2.2 x 10^^{M/10^'^MQf^'^-'^ergs s"^ . (3) 

The above relation is consistent with the theoretically expected dependence re- 
sulting from the thermal bremsstrahlung origin of the X-ray emission: La;(Bol) oc 
Mg^gT^-^ oc M^-^ (for an approximately constant size of the X-ray emitting re- 
gion, and Mgas oc M). This yields (0.1 — 2.4) oc M^, consistent with eq. 
(3). 

Using eq. (3) above, the observed X-ray flux of a cluster is 

(0.1 - 2.4keV) ^ 2.2 x 10^^ {M / 10^^ MQ^^^-^/And^ ergs cm'^ s"^ , (4) 

where d is the luminosity-distance of the cluster. The fliix is therefore propor- 
tional, approximately, to F oc (M/d)'^. For a flux-limited sample, the nearby 
clusters have a low mass (and thus a low richness) threshold, and the distant 
clusters have a richer threshold. In a flux-limited sample the number density 
of clusters decreases faster for poor clusters than for rich clusters. A volume- 
limited sample, on the other hand, has a cluster density that remains constant 
(in comoving coordinates) for each richness class. 

To calculate the expected correlation function of a flux-limited X-ray cluster 
sample, we use N-body simulations that match the correlation function of the 
observed richness threshold cluster samples (i? > 1 clusters, EDCC clusters, 
and APM clusters). We use a large-scale Particle-Mesh code with a box size 
of 400/i~^Mpc to simulate the evolution of the dark matter. The box contains 
500^ cells and 250^ = 10^-^ dark matter particles. The spatial resolution is 
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0.8/i~^Mpc. Details of the simulations are discussed in Cen (1992) and Bahcall 
& Cen (1992). A model that reproduces the observed mass-function of clusters, 
as well as the observed correlation function of the richness-threshold i? > 1, 
APM, and EDCC clusters is a low-density, unbiased CDM model (Bahcall & 
Cen 1992). This model, with = 0.2, h = 0.5, and no bias {b = 1) produces the 
observed richness-dependent universal cluster correlation function (eq. 1-2), and 
is consistent with other cluster properties such as their mass-function. 

Clusters are selected in the simulation box using an adaptive linkage al- 
gorithm following the procedure described in Suto, Cen & Ostriker (1992) and 
Bahcall & Cen (1992); the cluster mass within 1.5/i~-^Mpc is determined. In 
order to extend the volume-limited sample of the underlying cluster distribution 
to distances of the most distant X-ray clusters observed {z < 0.25), a mosaic of 
eight 400/i~^Mpc simulation boxes are used, corresponding to 800h~^ Mpc on a 
side. Within this larger mosaic box, ~ 3000 R > 1 clusters are identified in the 
simulation (with n = 6 x 10~^ h^Mpc~^). 

Each cluster of mass M (within 1.5/i~^Mpc) is assigned an X-ray luminosity 
as given by relation (3), and an X-ray flux — to an observer at the corner of the 
box — as given by (4). The flux threshold of the Romer et al. (1993) sample 
is then applied. All clusters in the 800/i~-'^Mpc simulation box with X-ray flux 
above this threshold are identified; they correspond to a fiiix-limited X-ray sample 
similar to the one observed. 

The redshift distribution of the X-ray clusters in both the observed and the 
simulated flux- limited samples are presented in Figure 1. Two simulated cases 
are shown: one corresponds to relation (4) (with a 15% lower amplitude in the 
Fx — relation in order to match the observed average co-moving density of 
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X-ray clusters in the range cz = 5000 to 50000 km s~-^, n ~ 5 x 10~^h^Mpc~^); 
the other corresponds to a somewhat shaUower slope than given in relation (4), 
Fx oc M^-"^, with an amplitude that yields an average comoving density of ~ 
8 X 10~^h^Mpc~^. Both simulations yield results that are consistent with the 
redshift distribution of the observed fliix-limited sample of clusters. Furthermore, 
we find that plausible variations in relation (4) do not produce significant changes 
in the results. 

The mean richness of the X-ray clusters as a function of redshift is presented 
in Figure (2). In the fliix-limited simulations, there is a strong increase of richness 
with redshift. By constrast, the mean richness in volume- limited samples is 
constant. Approximately 30% of the clusters are poor (i? < 0), consistent with 
the observed sample of Romer et al. (1993). 

The correlation function of the simulated flux-limited sample is presented 
in Figure 3. It is compared with the X-ray cluster observations of Romer et al. 
(1993). The agreement between the observations and simulations is excellent. 
The simulations yield f« 14/j,~^Mpc for the flux-limited sample, consistent with 
the observed To = 13.7 ± 2.3h~^Mpc. The correlation function of the simulated 
volume-limited rich R> 1 clusters (with n = 6 x 10~^h^Mpc~^, d = 55/t~^Mpc) 
is also shown, for comparison; this function matches well the observed R > 1 
correlations, with ~ 21/t~^Mpc (Bahcall & Soneira 1983; Postman et al. 1992; 
Peacock & West 1992). 

Figure 3 shows that the flux-limited sample exhibits a lower correlation am- 
plitude than the volume-limited, richness-limited sample of comparable average 
number density. This is expected due to the "richness-mixed" nature of the fliix- 
limited sample. In the present case, the correlation scales of the flux-limited and 
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the richness-limited samples satisfy 

2 

ro(X— ray fliix— limited) ^ -^foiR > 1) , (5) 

with ro{R > 1) ~ 21/i~-^Mpc and ro(X— ray fliix— limited) w 14/i~^Mpc, as 
observed. The results are insensitive to reasonable variations in the I'a; (-^cluster) 
relation. As an example, we present in Fig. 3 the correlation function for the 
case La,(0.1 - 2.4) = 3.5 x 10^2(M/10i^M©)i '^ergs s-\ which also matches the 
observed cluster density distribution (Fig. 1). 

3. CONCLUSIONS 

We show that the spatial correlation function of a flux-limited sample of X- 
ray selected clusters of galaxies exhibits a correlation scale that is smaller than 
the correlation scale of a volume-limited, richness-limited sample of comparable 
average number density (Fig. 3). This is expected due to the "richness- mix" 
nature of the flux-limited sample, which contains poor clusters nearby and rich 
clusters farther away. 

We show that the correlation scale of the new flux-limited X-ray cluster sam- 
ple from ROSAT (Romer et al. 1993) is expected to be ro(X— ray flux— limited) 
^foiR > 1) ~ 14/i~^ Mpc, as observed. We conclude that the new observations 
of X-ray clusters from ROSAT exhibit a correlation function that is consistent 
with, and is actually predicted from, the richness-dependent cluster correlation 
and ro{R > 1) 21h-^Mpc. 

It is a pleasure to acknowledge NCSA for allowing us to use their Convex- 
240 supercomputer, on which some of our computations were performed. R.Y.C 
is happy to acknowledge support from NASA grant NAGW-2448 and NSF grant 
AST91-08103 and N.A.B. acknowledges support from NSF grant AST93-15368. 
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FIGURE CAPTIONS 

Fig. 1- Redshift distribution of the fliix-limited X-ray clusters. The X-ray cluster 
observations of Romer et al. (1993) are shown by the faint histogram. Two 
simulated samples (§2) are represented by the dark and dashed histograms. 
dN{z) represents the relative number of clusters (i.e., fraction of total) in 
each redshift bin. (A volume- limited sample yields dN{z) oc z^). 

Fig. 2- The mean richness class of the flux-limited X-ray cluster samples as a func- 
tion of redshift. The two simulated samples (Fig. 1, §2) are shown. Also 
presented, for comparison, are the horizontal lines for the volume-limited, 
richness-limited R> 1 and R>0 clusters. 

Fig. 3- Cluster correlation function: the expected difference between the cluster 
correlation function of a volume-limited R > 1 cluster sample (faint line, 
with To ~ 21/i~-'^Mpc), and the relevant flirx-limited X-ray cluster correla- 
tions {Fx > 10~^^ergs cm~^ sec~^; dark and dashed lines, « 14/i~^Mpc). 
The flux-limited sample is expected to exhibit weaker correlations than the 
R > 1 complete sample. The observed flirx-limited X-ray cluster correla- 
tions (Romer et al. 1993), shown by the solid dots, are consistent with the 
expected flux-limited correlation function and ro{R > 1) ~ 21/j,~^Mpc. 
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